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^ Summary 

This  Report  considers  the  stress  analysis  of  elastic  plates  with  polar  orthotropic  properties.  The  analysis  is 
applicable  to  fibre  reinforced  composites  in  the  form  of  an  annulus  or  curved  strip  with  circumferential  fibres. 
For  certain  simple  loading  conditions,  including  those  of  pressure,  moment  and  shear,  curves  are  presented 
which  give  the  maximum  stresses  for  a wide  range  of  composite  geometries  and  elastic  properties. 


* Replaces  R.A.E.  Technical  Report  75060 — A.R.C.  36  239 
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1.  Introduction 

It  is  well  known  that  certain  fibres  have  exceptional  longitudinal  strength  and  stiffness  properties.  When 
such  fibres  are  embedded  in  a matrix  to  form  a 'fibre  reinforced  composite’  they  can  be  used  in  a structural 
context.  The  present  and  increasing  use  of  such  composites  is  paralleled  by  research  aimed  at  solving  the 
associated  structural  problems.  Considerable  progress  has  been  made,  for  example,  in  determining  the 
stiffness  and  stability  characteristics  of  rectangular  composite  plates'  ’,  and  in  determining  the  optimum  layout 
of  fibres  to  resist  specified  loads'  4. 

The  present  paper  focusses  attention  on  composites  with  curved  fibres.  Specifically,  the  analysis  relates  to  an 
elastic  plate  w ith  polar  orthotropy  such  that  under  conditions  of  plane  stress  or  plane  strain  there  are  constant, 
hut  generally  differing,  values  of  the  Young's  moduli  in  the  radial  and  circumferential  directions  E,  and  E„,  a 
constant  value  of  the  shear  modulus  G,„  and  constant  values  of  the  Poisson’s  ratios  i>,„  and  ver.  In  practice  this 
means  that  the  analysis  is  applicable  to  fibre  reinforced  composites  in  the  form  of  an  annulus  or  curved  strip 
w ith  circumferential  fibres  whose  density,  or  fibre  volume  fraction,  is  constant.  There  may  also  be  radial  fibres 
but  the  requirement  of  constant  density  means  that  such  fibres  should,  strictly  speaking,  be  of  varying  lengths 
to  accommodate  increasing  numbers  with  increasing  radial  distance.  Alternatively,  the  circumferential  fibres 
may  be  bonded  to  an  isotropic  plate  which,  itself,  may  be  composed  of  fibre  reinforced  laminates  in  a 
pseudo-isotropic  arrav  le.g.  0 degrees,  ±60  degrees).  The  analysis  does  not  cater  for  bending  of  the  plate,  so 
that  w here  such  bonding  occurs  the  arrangement  should  be  symmetrical  about  the  mid-thickness  plane.  Also, 
if  there  are  thickness  changes  in  the  plate  there  will  be  restrictions  arising  from  the  inherent  assumptions  of 
plane  stress  conditions,  furthermore,  such  restrictions  will  be  more  stringent5  than  those  in  the  corresponding 
isotropic  case  because  of  the  relatively  low  value  of  the  transverse  shear  modulus  Gri.  Finally,  in  the  context  of 
plane  strain,  the  analysis  is  applicable  to  a thick-walled  circular  cylinder  with  circumferential  fibre  reinforce- 
ment or  to  a unidirectional  fibre  reinforced  plate  formed  into  a sector  of  a cylinder,  e.g.  as  shown  in  Fig.  10.  In 
such  applications  the  modulus  E„  for  example,  refers  to  the  Young’s  modulus  through  the  thickness.  The 
multi-laminate  case  can  also  be  considered  provided  the  number  of  laminates  is  sufficiently  large  to  justify  the 
assumption  of  average  values  of  the  moduli  through  the  thickness. 

The  analysis  first  considers  the  simplest  stress  states,  namely  those  due  to  internal  or  external  pressure 
applied  to  an  annular  ring  or  thick-walled  circular  cylinder.  This  is  followed  by  the  pure  bending  of  a curved 
strip.  The  general  form  of  the  stress  function  is  next  derived  and  subsequently  used  for  the  analysis  of  the 
stresses  in  a curved  strip  subjected  to  a shear  force,  and  for  a problem  of  load  transfer  via  a semicircular  ring. 
All  but  this  last  problem  have,  indeed,  been  formally  solved  by  various  elasticians  of  the  Russian  school 
including,  in  particular,  I ekhnitskii  whose  book1  contains  a comprehensive  list  of  references.  The  motivation 
behind  the  present  (independent)  treatment  was  to  provide  an  engineering  insight  into  the  behaviour  of 
composites  with  polar  orthotropy.  The  treatment  accordingly  differs  from  earlier  work  in  the  extent  of  the 
numerical  results  presented.  Furthermore,  in  the  derivation  of  the  genera)  form  of  the  stress  function  certain 
differential  operators  are  introduced  which  markedly  simplify  the  subsequent  analysis. 
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For  more  general  types  of  loading  we  later  require  the  further  relation 


Now  if  r,„  is  zero  and  07,07,  are  independent  of  0,  as  in  the  present  problem,  the  equation  of  (radial) 
equilibrium  is 


07,-— (w,)  = 0, 

dr 


which  can  be  satisfied  by  the  introduction  of  a stress  function  'I7  such  that 


<r,~ — 

r 


By  the  same  token,  if  the  tangential  displacement  v is  zero,  as  in  the  present  problem,  the  direct  strains  can  be 
expressed  simply  in  terms  of  the  radial  displacement  u: 


The  equation  of  compatibility  is  therefore  given  by 


d 

er  = —(ree), 
dr 


which  can  be  expressed  in  terms  of  'F  in  virtue  of  equations  (1),  (2),  (5): 

2t  dV  ,d2V  n 
at  V-r—-r— r = 0, 
dr  dr" 


a2  = £,/£,. 

It  may  be  confirmed  that  equation  (8)  is  satisfied  by  terms  proportional  to  r±a,  and  it  follows  that  the  stresses 
can  be  expressed  in  the  form 

o’r  = —^T\  + Bp"\ 


a A . 

07,  = — —i+aBp  , 

P 


• 1 


where  p is  a nondimensional  measure  of  the  radius  given  by 


P = r/rl. 


(10) 
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2.1.  Internal  Pressure 

For  the  ease  of  uniform  internal  pressure  p i,  say,  the  constants  A,  H are  to  be  chosen  to  satisfy  the  conditions 


= ~P  I. 


and 


(11) 


WX-*  = o. 


where 


M ~r2/r 


(12) 


Thus 


A = 


(13) 


B = 


Of  particular  interest  is  the  maximum  value  of  the  hoop  stress  which  occurs  at  the  inner  boundary 

!>«]„. i d*) 

As/x  -»  1 the  hoop  stresses  increase  very  rapidly  and  in  Fig.  2 it  has  been  found  convenient  to  express  [cr*],,.,  as 
a multiple  of  the  nominal  hoop  stress  (r\P\/(r2-r\)),  i.e.  the  average  hoop  stress  required  to  equilibrate  the 
internal  pressure.  This  stress  ratio  can  therefore  be  regarded  as  a stress  concentration  factor  ai.d,  for  the  fibre 
reinforced  composite,  it  will  be  seen  that  for  given  values  of  r2/r , this  stress  concentration  factor  can  markedly 
exceed  that  for  an  isotropic  material. 


2.2.  External  Pressure 

The  case  of  uniform  external  pressure  p2,  say,  may  likewise  be  considered.  In  particular  the  hoop  stresses  at 
the  inner  and  outer  boundaries  are  given  by 


and 


/ Or  "*■  1 v 

k»]p-l  = -2orp:(  1,-  J, 


(15) 


In  the  isotropic  case  the  maximum  stress  occurs  at  the  inner  boundary,  but  for  the  fibre  reinforced  composite 
it  can  occur  at  either  boundary  depending  on  the  values  of  EJE,  and  r:/rt.  The  maximum  stress  has  been 
plotted  in  Fig.  3 as  a multiple  of  the  nominal  hoop  stress  (-r2p2/(r2- rt)).  When  the  maximum  occurs  at  the 
inner  boundary  the  corresponding  curve  is  shown  as  a broken  line:  when  at  the  outer  boundary  as  a full  line. 
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For  the  fibre  reinforced  composite  it  will  be  seen  that  for  given  values  of  r2/r , the  maximum  hoop  stress,  which 
generally  occurs  at  the  outer  boundary,  can  be  markedly  greater  than  that  in  an  isotropic  material. 
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2.3.  Equal  Internal  and  External  Pressure 

The  case  of  combined  internal  and  external  pressures  can,  of  course,  be  obtained  by  superposition  of  the 
results  of  Sections  2.1  and  2.2.  It  is  difficult  to  visualise  a practical  situation  where  this  might  occur  but  the 
results  for  the  case  of  equal  internal  and  external  pressure  p , say,  though  largely  of  academic  interest,  do 
highlight  the  differences  between  the  isotropic  and  fibre  reinforced  material.  In  the  isotropic  material  the  hoop 
and  radial  stresses  are  constant  throughout  and  equal  to  -p,  but  for  the  fibre  reinforced  material  (with  E„  * E,  I 
this  constancy  of  stresses  only  occurs  for  the  cut  ring.  For  the  continuous  ring  it  may  be  shown  that 


and 


= i =apy- 


=-«P 


-2m“+‘  + 1 
r“-l 


)■ 


( 16) 


' + 1 


M "I 


)J 


These  results  are  shown  in  Fig.  4 which  not  only  indicates  high  stress  concentrations  but  also  shows  that  the 
hoop  stress  at  the  inner  boundary  becomes  tensile  for  quite  modest  values  of  r2/ri  and  EJE,.  The  reason  for 
this  change  of  sign  and  the  high  stress  concentrations  in  this  and  the  previous  problems  is  readily  understood  by 
a consideration  of  the  limiting  case  in  which  E,  -*  0.  Then  the  applied  pressure  is  solely  resisted  by  the  adjacent 
circumferential  fibres  which  are  necessarily  very  highly  stressed. 


3.  Pure  Bending  of  Ring 

The  pure  bending  of  a cut  ring  or  curved  strip,  as  shown  in  Fig.  1 , is  also  characterised  by  zero  shear  stress  r,„. 
and  hoop  and  radial  stresses  which  are  independent  of  6.  But  the  tangential  displacement  v is  no  longer  zero 
and  the  equation  of  compatibility  must  be  derived  from  the  general  relationships  between  strain  and 
displacement.  Timoshenko*  shows  that  these  are 


e,  = = 


Hu 
Hr ' 


u 1 Hv 

ew=;+;^’ 


(a) 


(b) 


1 Hu  H (v\ 

Yr»  = ~ (-)  (0 

r rt 6 Hr  \r / 

In  the  present  problem  y,„  is  zero  and  hence,  from  equation  (17(c)), 

1 d2“_  a2  (v\ 

r2  HP2  HrHP\r) 

-Ur4 


(17) 


(181 


by  virtue  of  equation  ( 1 7(b)).  The  resulting  equation  from  which  v has  been  eliminated  may  now  be  combined 
with  equation  (17(a))  to  eliminate  u and  thus  derive  the  equation  of  compatibility: 


H e,  He,  H / 2Hee\ 

— t r — ^ — ( r — | = 0 
H6~  Hr  Hr\  Hr) 


(19) 


■ 


Substitution  of  equations  (1),  (2)  and  (5)  into  equation  (19)  enables  us  to  express  this  in  terms  of  the  stress 
function  '1/: 


a 


„ d2V  ,</>  n 
+ 2r— r + r'— T = 0. 


dr 


dr 


(20) 


It  may  be  confirmed  that  equation  (20)  is  satisfied  by  terms  proportional  to  r‘a  (as  before)  and  r.  Thus  the 
stresses  can  be  expressed  in  the  form 


and 


<ft 5TT  + Bp“  1 + C, 

P 


(21) 


aA  , 

o#  = — ZTi  + aBp  +C. 

P 


For  the  case  of  pure  bending  the  constants  A,  B,  C must  be  chosen  to  satisfy  the  following  conditions, 

l>r]P-i  =[ovW  =0,1 


and 

| raHdr  = - M, 

where  M is  the  moment  per  unit  length.  Hence  we  obtain 

MT[(n2a-Ha+i)  tt+1  „ 2„  . J 

<T’  = ~jr\ +(m  -1)p  — m +1  , 

ri  l p > 

and 

A/Tf  a(p2"-p“+1)  Q+I  2.  , .1 

2~l 0+1 + «(M  -l)p  ~P  +1  , 

r,  l p J 


(22) 


(23) 


where 


2(a2-  1) 

2a{(M2  + l)(p2“  + D-4M“  + '}-(a2+l)(p2-D(p2“-D' 


By  taking  the  limit  as  a -*  1 it  may  be  shown  that  these  results  agree  with  those  for  the  isotropic  case'’.  The 
maximum  tensile  and  compressive  hoop  stresses  occur  at  the  inner  and  outer  boundaries,  and  these  are  shown 
in  Fig.  5.  The  stresses  are  expressed  as  multiples  of  6 M/(r2-rl)2,  the  nominal  maximum  stress  assuming  a 
linear  hoop  stress  distribution.  This  form  of  presentation  enables  the  ordinate  in  Fig.  5 to  be  interpreted  as  a 
stress  concentration  factor.  It  will  be  seen  that  the  fibre  reinforced  material  exhibits  greater  stress  concentra- 
tions than  the  isotropic  material  but  the  increases  are  much  smaller  than  those  which  occur  in  a continuous  ring 
under  internal  or  external  pressure.  What  is  surprising  is  the  fact  that  while  the  maximum  hoop  stress  always 
occurs  at  the  inner  boundary  as  in  the  isotropic  case,  the  hoop  stress  at  the  outer  boundary  exceeds  the  nominal 
value  when  E„/E,  and  ri/ri  exceed  certain  critical  values.  The  reason  for  this  becomes  clearer  if  we  look  at  Fig. 
6 which  shows  the  radial  distribution  of  hoop  stress  in  rings  in  which  r2lr,  = 2 and  Ee/Er  = 1 or  50.  (Results  for 
intermediate  values  of  E„/E,  are  adequately  given  by  linear  interpolation.)  Compared  to  the  hoop  stresses  in 
the  isotropic  ring  those  in  the  fibre  reinforced  ring  exhibit  a superimposed  S-ing  distribution  which  accentuates 
the  stresses  near  the  inner  and  outer  boundaries  at  the  expense  of  those  towards  the  mid-thickness. 
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The  radial  stress  reaches  a maximum  value  at  a radius  specified  by 


r 


* 

I 

i*: 


(24) 


whence,  from  equation  (23), 

p=p*,  say 

The  maximum  radial  stress  a*,  say,  given  by  equations  (23)  and  (25),  is  shown  in  Fig.  7 as  a multiple  of 
6A//(r2  — 4 - 1)/ (/u. ^ + 1)].  The  introduction  of  the  additional  factor  [(/a*  - l)/(p*  + 1)]  means  that  cr*  is 

expressed  as  a multiple  of  the  nominal  maximum  assuming  a linear  hoop  stress  distribution.  [Briefly,  if 


1 +/a  -2p\ 

M-l  r 


it  follows  from  equation  (4)  that 


V = ~ \P  cre  dp, 
P h 

= ^l(M 


+ 1 ~P 


whose  maximum  value  is  [(/x 4 - 1)/(m*  + 1)].] 

It  will  be  seen  that  the  maximum  radial  stress  is  somewhat  less  in  the  fibre  reinforced  case,  though  the 
reduction  is  very  small  when  r2/r\  < T5. 


3.1.  Fibre  Failure  versus  Interlaminar  Failure 

In  a ring  under  pure  moment  the  maximum  value  of  the  radial  stress  is  much  smaller  than  the  maximum 
hoop  stress,  and  in  an  isotropic  material  the  radial  stress  has  no  bearing  on  failure  or  the  onset  of  plasticity. 
However,  in  a fibre  reinforced  material  with  circumferential  reinforcement  alone  the  radial  stress  necessary  to 
cause  failure  of  the  matrix  or  bond  is  much  smaller  than  the  hoop  stress  necessary  to  cause  failure  of  the  fibres. 
Thus  it  is  of  interest  to  determine,  for  example,  the  critical  dividing  line  between  such  modes  of  failure. 

The  ratio  of  the  maximum  hoop  stress  [cr«]p,i,  denoted  by  cr*,  to  the  maximum  radial  stress  a*  is  shown  in 
Fig.  8 for  various  values  of  r2/r,  and  EJE,.  For  unidirectional  CFRP,  for  example,  typical  values  for  the  ratio 
of  the  failing  stresses  under  longitudinal  affd  transverse  tension  lie  in  the  range  20-30.  If  we  confine  attention 
to  values  of  r2/r,  which  yield  values  of  <t*/ct*  in  this  range  it  is  seen  that  the  influence  of  the  modulus  ratio 
E„/E,  is  very  small.  In  particular,  if  the  ratio  of  failing  stresses  is  20,  say,  it  follows  from  Fig.  8 that  interlaminar 
tensile  failure  will  occur  before  tensile  fibre  failure  if  r2/rx>\- 25,  while  if  the  ratio  is  30,  say.  the 
corresponding  range  is  r2/rx  > T15. 
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3.2.  The  Flexural  Rigidity  of  the  Ring 

We  conclude  our  discussion  of  the  pure  bending  of  a cut  ring  with  a determination  of  the  effective  flexural 
rigidity.  In  this  connection  it  is  convenient  to  express  this  as  a multiple  of  the  flexural  rigidity  of  a straight  beam 
of  modulus  Efl,  depth  (r2  — r,)  and  unit  width,  the  length  of  this  straight  beam  being  i(r,  + r2)0. 

The  simplest  way  to  determine  the  effective  flexural  rigidity  is  via  the  displacement  v,  which  is  the  only 
displacement  component  which  contributes  to  a rotation  of  sections.  By  the  same  token  it  follows  from  Section 
2 that  the  A and  B components  of  the  stresses  in  equation  (2 1 ) do  not  contribute  to  the  rotation  and  we  are  left, 
somewhat  surprisingly,  with  the  constant  stress  component  C.  The  ensuing  analysis  is  accordingly  simpler  than 
that  for  the  isotropic  case! 
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If  the  index  c is  introduced  to  identify  this  constant  stress  component  we  can  write 


C C /O 

cr,  = a„=  C, 


(26) 


and  equation  (1)  shows  that  the  corresponding  strain  components  ec,  and  F»are  constants,  and  hence  equations 
(17)  may  he  readily  integrated  to  yield 

u = rO^e#— + rigid  body  terms.  (27) 

Now  the  rotation  it  of  a section  at  0 is  given  by 


(28) 


and  hence  the  change  in  curvature  of  the  mid-thickness  line,  whose  length  is  specified  by  s = 2<r,  +r2)8,  is  given 
by 


dn 

ds 


2(e;-ec„) 

r,+r2 

2C(a2—  1) 
Ee(r,+r2)  ’ 


(29) 


in  virtue  of  equations  (1)  and  (26).  Also,  in  terms  of  an  effective  flexural  rigidity  J,  say,  engineers’  theory  of 
bending  gives 


dn_M 

ds  ~ J’ 


(30) 


and  hence  from  equations  (23),  (29), 


J _ 6(1  +/*) 

£^~r(M  — i)V:-i)(i— m;") 

where 


(31) 


The  variation  of  J/EJ  with  r2/r,  for  various  values  of  EJE,  is  shown  in  Fig.  9.  It  is  seen  that  for  the  isotropic 
case  the  effective  flexural  rigidity  is  adequately  given  by  El  for  values  of  r2/r , up  to  about  3.  In  contrast,  the 
fibre  reinforced  case  shows  a marked  reduction  in  J for  all  but  the  smallest  values  of  r2/r,. 


4.  General  Loading  of  Ring 

The  loadings  considered  so  far  have  been  characterised  by  the  vanishing  of  one,  or  both,  of  the  terms  tv#  and 
v,  and  this  has  resulted  in  simplified  equations  of  equilibrium  and  compatibility.  In  the  general  case  the 
equations  of  equ'Hbrium  are6 


d<rr  j 1 dTrg  ^ (Tr-(T„ 

dr  r dO  r 


and 


(32) 


1 da-„  | ^ 2rr9 

r dO  dr  r 
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and  these  are  satisfied  hy  the  introduction  of  a stress  function  <1>  such  that 


I 


(Jr 


1 M>  1 r’<l> 

+“ S TT  t 

r dr  r <W~ 


<rH 


Trv 


a1® 

dr 2 

a / l ^ 

dr\r  dd 


(33) 


As  for  the  equation  of  compatibility,  it  may  be  confirmed  by  substitution  that  the  general  strain  displacement 
relations  of  equation  (17)  satisfy  the  equation 


d 

1 

1 dy,9\) 

dr 

l r V ar " 

" r ae  ) 1 

d2er 

~d(P 


0. 


(34) 
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f 


**» 

j 


» 

1 

C* 

f 

* 


i 

r* 

* 
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The  equation  satisfied  by  the  stress  function  $ is  obtained  by  substitution  of  equations  (1),  (3)  and  (33)  into 
equation  (34).  This  stress  function  equation  is  best  expressed  in  a compact  form  by  introducing  the  following 
differential  operators  defined  by: 


and 


(35) 


It  may  now  be  shown  that 

[*  A?  + 2/8  A,2  A,  + M - ( 1 + k + 2/3 ) A;]<D  = 0,  ] 


where 


k = E'/E, 


and 


(8  = 


E, 

2 Gre 


«v 


Note  that  for  the  isotropic  case 


k =0  = 1, 

and  the  differential  operator  in  equation  (36)  factorises  to  give 

[Ar2  + A8±2Ar]<t>  = 0, 


(36) 


(37) 


whence 


= +<J>2,  say 
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where 


f 


V’<t>,=0, 


and 


1 , ; 

<I>2  - r— - + - r V <l>2  = 0, 
dr  4 


(38) 


and  we  recover  the  biharmonic  equation,  albeit  in  a separated  and  unusual  form. 

4.1.  General  Form  of  the  Stress  Function 

A general  expression  for  <t>  can  readily  be  obtained  in  the  form  of  a series 


4>=  I 

n =0 


(39) 


where 


41,,  = X (a„,r'4‘  cos  nO  + bn,rq‘  sin  nd), 


(40) 


where  the  ani  and  bm  are  constants  and  the  indices  q,  are  determined  below.  Substitution  of  equation  (40)  into 
equations  (35)  and  (36)  shows  that 


A,<t>„  = (<? -!)<!>„ 


and 


=(l-rr)<J>„, 


(41) 


and  it  follows  that  for  any  given  value  of  n there  are,  in  general,  four  values  of  q (denoted  by  q, ) determined  by 
the  equation 


/c(<7-1)4-(4-1)j(1  + k + 2/?«2)  + («2-1)2  = 0. 


(4: 


When  either  q or  n equals  unity  equation  (42)  results  in  repeated  roots  and  further  analysis  is  required.  Thus  if 
q = 1 we  must  consider  the  general  expression 


<t>  = rf(0). 

Substitution  into  equations  (35)  and  (36)  shows  that 

Arct>  = 0, 


(43) 


so  that 


a;/  = o, 


which  may  be  integrated  to  give 


f(0)  - a\8  sm  6 +b]0  cos  0 + (terms  proportional  to  sin  8 and  cos  0.  which  do  not  contribute  to  the  stresses). 

(44) 


By  the  same  token,  if  n = 1 we  consider 


4>.=«(r)(*,ny 

Vcos/ 


(45) 
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...  \ 


Substitution  into  equation  (35)  shows  that 


A»<1>|  = 0, 


and  hence  equation  (36)  may  be  expressed  in  the  form 

£(r)  = gi(r)  + g2(r),  (46) 


where 


[*Af2 * * * *-(l+*  + 2,3)]g,=0, 


(47) 


and 


tigz  = 0.  (48) 

Equation  (47)  admits  of  a solution  consistent  with  equation  <4d),  while  equation  (48)  may  be  integrated  to  give 

) = cr  In  r + (a  term  proportional  to  r which  does  not  contribute  to  the  stresses).  (49) 

Finally,  consideration  of  the  limiting  case  in  which  n = 0 leads  us  to  the  solution 

d>„  = a o.  + box  e + (a02  + bU2e)r ,+“  + (ani  + bm8)r'a  + (aU4  + bU48)r\  (50) 

where  a is  as  defined  in  equation  (8).  In  this  expression  for  <b<,  the  term  aU\  does  not  contribute  to  the  stresses, 
while  the  terms  identified  by  a02,  um  and  correspond  to  the  stress  function  'F  of  Sections  2,  3. 

[Note  that  each  of  the  solutions  derived  in  this  section  has  its  counterpart  in  the  general  stress  function  for 
the  isotropic  case  which  was  first  derived  by  Michell7  and  is  reproduced  in  Ref.  6.  The  present  analysis, 
however,  shows  that  two  solutions  were  overlooked  in  this  earlier  work,  namely 

<t>  = A0  \nr  + Bdr2  In  r.] 


2.  Shear  Force  applied  to  Curved  Strip 

The  pure  bending  of  a cut  ring  or  curved  strip  was  considered  in  Section  3.  This  analysis  is  augmented  here 

by  a consideration  of  the  stresses  in  a curved  strip  subjected  to  shear.  A situation  in  which  this  would  be 
relevant  is  shown  in  Fig.  10  which  depicts  a section  through  the  edge  region  of  a fibre  reinforced  plate  A where 
it  is  attached  to  another  member  B at  right  angles.  The  edge  of  the  plate  has  been  formed  to  include  an  annular 
quadrant  and  any  transverse  loads/pressures  acting  on  the  plate  may  be  resolved  at  section  C into  a moment  M 
per  unit  length  and  a shear  force  S per  unit  length.  At  section  D the  shear  force  S is  resisted  by  a moment  and  a 

force  S which  is  now  parallel  to  the  fibres;  the  solutions  for  a moment  and  shear  force  can  therefore  be 

combined  to  yield  also  the  solution  for  a force  applied  in  any  direction. 

The  stresses  due  to  the  shear  force  S may  be  derived  from  the  stress  function  of  equation  (45)  which,  for 
con'  lienee,  we  express  in  the  form 


where 


and 


cp,  = sin  *|  + (^y  + a2p  In  P j, 

<?1.2=  1±0>. 


(51) 
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The  derived  stresses  are  given  by 


* 

♦ 

\ 


;y 


1 

r- 


i 

- 


cr,  = p 2 sin  0(a]p‘l'  + a2p'l2  + a^p), 
cr„=p  2 sin  0(q,a  1p,<l  + q2a2p'*-'  + atp), 
and 

t,»  = -p  2 cos  0(a ip*1  + a2p‘l:  + d)p). 


and  it  is  to  be  noted  that 


(52) 


ov/sin  0 = -T,.<,/cos  6. 

The  boundary  conditions  along  the  curved  edges  are  that 

frrr  Trtt  0 Jr  — r i ,r_.  * 

while  equilibrium  at  section  C requires  that 


-r 


T*dr. 


Equations  (54)  and  (55)  suffice  to  determine  the  constants  a,  which  are  given  by 


S Sp " S(p“  + 1) 

a\-  „ - ai~  r^andflj—  — , 

r\K  r{K  r^K 


where 


2(a"  — 1) 

K = — ——(p“  + 1)  In  p. 


(53) 


(54) 


(55) 


(56) 


It  is  seen  from  equation  (56)  that  the  elastic  properties  affect  the  stresses  only  insofar  as  they  affect  the 
parameter  w.  More  directly  it  is  seen  that  the  governing  parameter  is  given  by 


1+2/3 
r,,say  = — 


lh 

r (±  + ±\ 

3 i 

(57) 


where  the  factor  1/3  has  been  introduced  so  that  r/  = 1 in  the  isotropic  case. 

The  variation  of  the  hoop  and  shear  stresses  through  the  thickness  of  the  curved  sheet  (or  beam)  is  shown  in 
Figs.  1 1 , 1 2 for  the  case  in  which  r2/rx  = 2 and  r/  = 1 or  50.  If  the  hoop  stresses  are  compared  with  those  in  Fig. 
6,  which  relates  to  an  applied  moment,  it  is  seen  that  for  equal  values  of  EJE,  and  rj  the  hoop  stresses  in  the 
fibre  reinforced  case  exhibit  a greater  degree  of  S-ing  due  to  an  applied  shear  force.  This  is  reflected  in  Fig.  12 
in  a flatter  distribution  of  the  shear  stress. 

The  maximum  value  of  the  hoop  stress  occurs  at  0 = {ir,  r = ru  and  is  given  by 


&9.max 


Sa»(l-p“) 

rtK 


(58) 


which  is  shown  in  Fig.  13  where  oyma„  has  been  expressed  as  a multiple  of  r2S/(r2-r])2  in  order  to  flatten  the 
curve  for  the  isotropic  case  and  so  high-light  the  influence  of  anisotropy. 


13 


j 


The  shear  stress  and  the  radial  stress  reach  maximum  values  at  a radius  specified  by 


r 


p**,  say 


/ 2Uo  + l)  \1/" 

' “ + 1 )"  + 4a  (w~  — l ))*  + 11“  + 1 / 


(59) 


and  Fig.  14  shows  the  corresponding  values  of  TrU,ma,  expressed  as  a multiple  of  -S/(r2  - ri).  It  is  seen  that  as 
r2/r i -»  1 this  stress  ratio  approaches  T5  which  is  the  value  appropriate  to  a parabolic  distribution  of  shear. 
Deviations  from  this  value  are  comparatively  modest.  Fig.  14  is,  of  course,  also  applicable  to  <j,  nlJX  in  virtue  of 
equation  (53). 


4.3.  Load  Transfer  via  Semi-circular  Ring 


As  a further  example  in  the  use  of  the  general  stress  function  we  consider  a purely  circumferentially  fibre 
reinforced  semi-circular  ring  with  contiguous  fibre-reinforced  parallel  strips  as  shown  in  Fig.  15a.  A force  P 
per  unit  length  is  applied  along  the  line  0 = 0,  via  a pin  of  radius  r,,  and  is  resisted  by  loads  of  magnitude  \P  in 
each  of  the  strips.  Strictly  speaking  this  is  a difficult  contact  problem  which  requires  a knowledge  of  the 
clearance  and  frictional  properties  between  pin  and  ring,  and  the  elasticity  and  support  conditions  of  the  pin. 
etc.  As  such  it  is  outside  the  scope  of  the  present  paper  but  the  essential  features  may,  nevertheless,  be 
reproduced.  Thus,  to  simplify  the  problems  we  will  assume  that  along  the  interface  between  pin  and  ring  the 
radial  stress  varies  as  cos2  0,  so  that 


3 P 


l>r]r»r,  = “jjHl  + COS  20), 
or  i 


and 


[>,»],=,,  = 0, 


(60) 


while  across  the  sections  at  8 = ±\ir. 


and 


le  = *i 


. = 0. 


(61) 


These  boundary  conditions  are  the  same  as  those  in  a complete  ring  subjected  to  ‘mirror  image"  loads  as  shown 
in  Fig.  15b.  The  solution  corresponding  to  the  constant  stress  component  in  equation  (60)  is  given  in  Section 
2.1  while  the  solution  for  the  part  which  varies  as  cos  20  can  be  obtained  from  the  stress  function 


= cos  20  £ apq\ 

i-1 


(62) 


where  the  indices  q,  are  the  roots  of  the  equation 


K(q-  1)4-(<?-1)2(1+k+8/?)  + 9 = 0, 


(63) 


and,  for  convenience,  we  have  reintroduced  the  nondimensional  term  p(=  r/r\). 
Substitution  of  equation  (62)  into  equation  (33)  gives 


<r,  = cos  20  £ (q,  - 4 )ap 


4 ,-2 


ov,  = cos  20  £<f,(<f,  - Dap"1  2, 


and 


Tre  = sin  20  I (q,  -\)apq‘  2. 


(64) 
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The  boundary  conditions  of  equation  (61)  are  automatically  satisfied  by  the  expression  for  <t>,,  while  the 
vanishing  of  07  and  t,„  at  r = r:  (i.e.  p = p.)  implies  that 


L • 

j. 

L 


Za,p*-  = 0, 

and  • ^5) 

Similarly  equations  (60)  can  be  combined  to  yield 

?*-£• 

and  (66) 

£(<?,-  l)a,  =0. 

/ 

Equations  (65)  and  (66)  suffice  to  determine  the  coefficients  a,  and  hence  the  maximum  value  of  o-„,  which 
occurs  at  the  points  r = ru  0 = can  be  determined  from  the  relation 

o'*,  say  = I «?,(<?,-  l)a,.  ,67) 

As  for  the  constant  stress  component  in  equation  (60)  it  will  be  seen  by  comparison  with  equations  ( 1 1 1 and 
(14)  that  the  maximum  value  of  cr«  is  given  by 


trj*,  say  = 


7>aP  I p2a  + 1( 


8 r,  \p 


Thus,  in  the  problem  under  consideration  the  maximum  value  of  is  given  by 

(To.  max  = 0-J+O-J*.  (69) 

Calculations  have  been  made  for  the  following  elastic  parameters: 

0 = 0-5,  1 and  k — 1/30. 

The  values  chosen  for  0 are  effectively  lower  and  upper  bounds  for  all  unidirectional  fibre  reinforced 
composites,  while  the  value  chosen  for  k is  typical  of  CFRP.  Table  1 below  shows  for  various  values  of  r2/ru 
the  corresponding  values  of  ov,  max  expressed  as  a multiple  of  the  nominal  hoop  stress  P/2(r2  — r ,).  These  stress 
ratios  can  therefore  be  regarded  as  stress  concentration  factors. 


TABLE  1 

Valuesof^.m.«/(^r^j) 


r2/r,  = 1-25 

1-5 

2-0 

0 =0-5,»c  = 1/30 

9-69 

8-51 

10-83 

0 = 1,  k = 1/30 

10-34 

9-51 

12-07 

Before  discussing  these  results  it  is  important  to  remember  that  they  apply  strictly  to  the  load  system  depicted 
in  Fig.  15b.  Some  general  conclusions  will  now  be  drawn:  The  stress  concentration  factors  are  high  but  vary 
with  r:/r,  by  relatively  small  amounts;  this  is  because  as  r2/r1  increases  the  component  of  the  stress 
concentration  factor  due  to  the  constant  pressure  term  cr* * increases,  while  that  due  to  the  bending  component 
<rj  decreases.  Also,  as  0 varies  from  0-5  to  1 there  is  an  increase  in  the  stress  concentration  factor  of  about  10 
per  cent.  Some  such  increase  is  to  be  expected  because  this  change  in  the  value  of  0,  with  k constant,  implies  a 
reduction  of  about  50  per  cent  in  the  shear  modulus  Gr„. 

Note  that  recent  work  suggests  that  the  pressure  distribution  between  the  pin  and  ring  is  approximately 
constant.  If  this  is  so,  it  can  be  shown  that  the  stress  concentration  factor  increases  monotonically  with  r2/rt, 
and  at  r2/r{  = 2-0,  k = 1/30,  for  example,  it  is  equal  to  about  4-0. 

5.  Conclusions 

This  Report  considers  the  stress  analysis  of  elastic  plates  with  polar  orthotropic  properties.  It  is  applicable  to 
fibre  reinforced  composites  in  the  form  of  an  annulus  or  curved  strip  with  circumferential  fibres;  in  some 
circumstances  it  is  also  applicable  to  a multi-laminate  composite.  The  simple  loading  conditions  of  pressure, 
moment  and  shear  have  been  considered  in  detail  and  curves  are  presented  which  give  the  maximum  stresses 
and  other  features  for  a wide  range  of  composite  geometries  and  elastic  properties.  As  a general  rule  the  stress 
concentration  factors  in  the  composite  are  higher  than  those  in  the  corresponding  isotropic  case,  though  there 
are  minor  exceptions.  In  the  derivation  of  the  general  form  of  the  stress  function  certain  differential  operators 
are  introduced  which  markedly  simplify  the  subsequent  analysis. 
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/ (r2-ri)Vl2 

J Effective  flexural  rigidity 

M Moment  per  unit  length 

p Pressure 

P Applied  force  per  unit  length 

q Index 

r , 6 Polar  coordinates 

r i,  r2  Inner  and  outer  radii 

S Shear  force  per  unit  length 

u,  v Radial  and  tangential  displacements 

a,  f),  k,  w,  q Parameters  defined  by  equations  (8),  (36),  (5 1 ),  (57) 

X , X Differential  operators  defined  by  equation  (35) 

e„  £».  yre  Direct  and  shear  strains 

'I',  <t>  Stress  functions  defined  by  equations  (5),  (33) 

M 'i/'i 

i/r„,  v„r  Poisson’s  ratios 

P 'Ai 

cr„  <r„,  t,„  Direct  and  shear  stresses 

fl  Rotation  of  section. 
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Fig.  11.  Variation  of  hoop  stress  through  thickness  (rs/r,  = 2). 
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